Field-induced gap in ordered Heisenberg antiferromagnets 
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Heisenberg antiferromagnets in a strong uniform magnetic field H are expected to exhibit a gap- 
less phase with a global 0(2) symmetry. In many real magnets, a small energy gap is induced by 
additional interactions that can be viewed as a staggered transverse magnetic field h = cH, where 
c is a small proportionality constant. We study the effects of such a perturbation, particularly for 
magnets with long-range order, by using several complimentary approaches: numerical diagonaliza- 
tions of a model with long-range interactions, classical equations of motion, and scaling arguments. 
In an ordered state at zero temperature, the energy gap at first grows as (cH)^^^ and then may dip 
to a smaller value, of order (cH)^^^ , at the quantum critical point separating the "gapless" phase 
from the gapped state with saturated magnetization. In one spatial dimension, the latter exponent 
changes to 4/5. 



I. INTRODUCTION 

The investigation of quantum antiferromagnets in 
strong magnetic fields is currently a very active field of 
research. Several remarkable properties have been ob- 
served, such as the closing of the energy gap in spin-1 
chains^ and spin ladders^ and observations of magneti- 
zation plateaux in frustrated magnets. '^•^ While the broad 
features of these models can be explained in the frame- 
work of the Heisenberg model in a uniform magnetic field, 
a closer examination of experimental data reveals devia- 
tions from theoretically predicted behavior. In particu- 
lar, the supposedly gapless phases actually possess small 
energy gaps, which can only be explained by the presence 
of anisotropic interactions. 

For instance, Dender et al^ discovered that an applied 
magnetic field H induces in spin- 1/2 Heisenberg chains 
a gap A c>c 770-65(3). Oshikawa and Affleck^ ascribed the 
gap to a staggered transverse field arising from the stag- 
gering of the g-tensor or of the Dzyaloshinskii-Moriya 
(DM) interaction. They suggested the effective Hamilto- 
nian 



Ti. = ^[JS„ • S„+i — HS^ 



(1) 



where h <x. H . The transverse field h creates a spin 
gap A cx h?^^ (X H"^/^ , in agreement with the experi- 
ment and density-matrix renormalization group (DMRG) 
calculations 

An extension of these results to higher dimension is 
still at a preliminary stage. Sato and OshikawE^ studied a 
model of interacting chains with a staggered field induced 
by the DM interactions and found that the gap scales 
as A (X H^/^ in a weak field. The recent discoverji^ 
of the Shastry-Sutherland antiferromagnet SrCu2 (603)2 
calls for further studies in two dimensions. In particular, 
the presence of a staggered magnetization in the low- 
field NMR signal was ascribed to a staggering of the g- 



tensor and of the DM interaction.-^ The direct numerical 
investigation of the relevant microscopic models does not 
seem to be possible at the moment. Indeed, to get a 
reliable estimate of the small gap induced by anisotropic 
interactions requires to reach sizes such that the finite- 
size gaps, typically of order J/N for N sites, are much 
smaller than the physical gap. This is possible in ID with 
the help of the DMRG algorithm, which by now routinely 
allows one to study systems with 200 sites or more, but 
not in 2D. 



In this paper, we present the first systematic study of 
this problem in the context of the Lieb-Mattis modeljii 
wherein every spin of one sublattice is coupled equally 
to all spins of the other. This model is expected to pro- 
vide a fair description of the long-wavelength properties 
of bipartite Heisenberg antiferromagnets with long-range 
Neel order in the ground state. Therefore it should be 
relevant for magnets in 2 and 3 dimensions. 



We pay particular attention to the high-field regime 
where the uniform magnetization saturates. In the ab- 
sence of a staggered field, the saturation occurs at a crit- 
ical point H — He separating a gapless ordered phase 
with a spontaneously broken rotational 0(2) symmetry 
from a gapped phase with fully polarized spins. The 
behavior of the transverse-field spin gap near the criti- 
cal field He is an important problem in view of its rele- 
vance to a number of experimental systems. Unlike the 
weak-field regime, for which analytical results have been 
obtained, the saturation region has so far been studied 
only numerically.'''^ Curiously, the numerical data reveal 
a nonmonotonic dependence A(iJ), with a pronounced 
minimum near the saturation field He- A similar effect 
was noticed earlier by Sakai and Shiba in their numerical 
analysis of spin-1 chains 
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II. THE LIEB-MATTIS MODEL 

In this section, we concentrate on the Lieb-Mattis 
model describing a Heisenberg antiferromagnet on a bi- 
partite lattice in which every spin of one sublattice is 
coupled through antiferromagnetic exchange to all sites 
of the other sublattice. Its Hamiltonian describes spins 
of length S" in a uniform magnetic field H = (0, 0, H) and 
a staggered field h = {h, 0, 0) perpendicular to it: 



H = [J/N) Y^Y.^-- 
rSAr'eB 

-5](H + h).Sr-^(H-h).Sr (2) 

= {J/N) Sa-Sb- H{S\ + 5|) - h{Sl ~ SI), 

where Sa and S b are the total spins of sublattices A and 
B. The exchange constant J is normalized by the total 
number of sites N to make the energy of the model an 
extensive quantity 0{N). To reflect the induced origin 
of the staggered field h, we will set h = cH, where the 
constant c ^ 1 is a property of the system. 

In zero transverse field, the Hamiltonian ^ has an 
0(2) rotational symmetry and is readily diagonalized. 
The energy levels are expressed in terms of the total spins 
of the sublattices Sa and Sb , the total spin S'tot and its 
projection S^^^ on the direction of the uniform field H: 

Eo = {J/2N)StotiStot + l)-HS^,^ 

- iJ/2N)[SA{SA + l) + SB{SB + l)]. (3) 

In the ground state the sublattices are fully polarized, 
Sa — Sb — NS/2, and the total angular momentum 
points along H: S^^^ — Stot- Excitations reducing sub- 
lattice magnetizations have an energy gap 0{JS). Below 
saturation, H < He = JS, low-energy excited states are 
obtained by changing the quantum numbers 5*^^^ = S'tot 
from their ground-state values. These excitations have 
energy 0{JS/N) and form a continuum in the limit 
N oo. The gapless excitations are caused by spon- 
taneous breaking of the 0(2) symmetry. Indeed, for 
N —^ oo the sublattice spins and become classi- 
cal vectors with well-defined directions in space (Fig.^). 
In the absence of the transverse staggered field h, the 
sublattice moments can be freely rotated about H. The 
low-energy excitations thus correspond to a slow preces- 
sion of Sa and Sb when the angle 9 deviates slightly 
from the equilibrium value = arcsin {H/ JS) . 

Adding the transverse field h breaks the 0(2) symme- 
try explicitly and violates conservation of both 5tot and 
Sjoj. Nonetheless, lengths of the sublattice spins Sa and 
Sb are still good quantum numbers. This fact greatly 
simplifies the numerical diagonalization of the Hamilto- 
nian 0. 
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FIG. 1: A classical ground state. The sublattice spins reside 
in the plane formed by the vectors H and h. 

A. Exact diagonalizations 

Because the ground state and all low-lying excitations 
of the model with h = are in the sector Sa — Sb — 
NS/2, we will restrict our analysis of the gap to that 
sector. The size of this subspace for N spins equals 
{NS -1-1)^, which is smaller than the size of the total 
Hilbert space, (25* -1-1)^. This enables us to treat systems 
with rather large numbers of spins. In the following, we 
present results for 2000 spins S = 1/2, a size clearly be- 
yond the scope of exact diagonalizations of other Heisen- 
berg models without additional conserved quantities. 

We first calculate the gap as a function of h in the 
absence of a uniform field. The results are plotted in 
Fig|21 By fitting the data at low field h ^ 0, we deter- 
mined that the gap vanishes as A ~ vJSh, precisely as 
found by Oshikawa and Affleck in the approximation of 
noninteracting magnonsi^ This result is most easily un- 
derstood by computing the precession frequency of sub- 
lattice magnetizations at the classical level, as we do in 
the next section. The result of this calculation is plotted 
as a solid line in Fig. [21 An excellent agreement shows 
that size effects are already negligible for N = 2000 (with 
the exception of a finite gap at h = 0). 

Next we consider the case where the staggered field 
is proportional to the uniform field: h = cH. Typi- 
cal results obtained for various values of c are plotted 
in Fig. 121 As 0, A ~ H^^^ with numerical preci- 

sion. The behaviour at higher fields depends on the value 
of the proportionality constant c. If it is small enough, 
c < 0.03 ± 0.005, the gap exhibits a local minimum close 
to — and slightly below — the saturation field. Such a dip 
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FIG. 2: The spin gap A as a function of the transverse field 
h for two values of the uniform field: H = (triangles) and 
H = He = JS (squares) for N = 2000. The curves are the 
results of the spin-wave analysis. 



FIG. 4: Uniform (filled squares) and staggered (open squares) 
magnetizations for A'' — 2000 as a function of H for c — 0.03. 
Solid and dashed lines are results of the spin-wave analysis. 
Dash-dotted lines show the behavior for h = 0. 



0,4 




H/JS 

FIG. 3: The spin gap as a function of at ft = cH for several 
value of the proportionality constant: c = 0.1 (squares), 0.03 
(triangles up), 0.01 (triangles down) and 0.0025 (circles). The 
gap has a local minimum for c < 0.03. Lines are the results 
of the spin-wave analysis. 



has been observed previously in numerical studies of ID 
models;^ but no explanation has been given so far. To 
get further insight, we have kept the uniform field at its 
saturation value He = JS and calculated the gap as a 
function of h. The results are plotted in Fig [3 They are 
consistent with a power-law scaling A cx h^^^. We will 
return to this phenomenon in the next section. 

Finally, we have calculated the uniform and staggered 
magnetizations as functions of H for several values of 
c (see Fig. 0}. As expected, the introduction of the 
symmetry-breaking staggered field removes the critical 
behavior near H = He = JS. 



B. Semi-classical analysis 

The results of the exact diagonalizations can be under- 
stood in the framework of the classical equations of mo- 
tion for the sublattice moments Sa and Sb. The frequen- 
cies of classical spin waves are identical to the energies 
of magnons obtained in the lowest-order 1/S expansion, 
one of the methods employed by Oshikawa et al&^'i^ Be- 
cause the magnon gap is lowest at zero wavevector, we 
specialize to uniform solutions, and consider precession 
of sublattice spins Sa and at maximal length NS/2. 

The equilibrium canting angle 6 (Fig. ^ is determined 
by minimization of the classical energy: 

JS'^ sine cose + hS sine -HS cose = (4) 

In the absence of the transverse field, 9 = arcsin (H/Hc) 
for H < Hc^ JS and 7r/2 for H > He- For h ^ 0, the 
kink in 0{H) is smoothed out. 

The equations of motion for the sublattice moments 
Sa and Sb arc: 

TlSa = (JSs/iV-H-h) X Sa, 

TiSb = {JSa/N -U + h) xSb- (5) 

We rewrite these equations in terms of the uniform and 
staggered magnetizations m — {Sa + Sb)/N and n = 
(S^ — Sb)/N and linearize them in small deviations Sm 
and Sn from the equilibrium values: 

hSih = — H X 5m — h X (5n, 
nSh = -(Jn + h) X (5m + (Jm - H) X 5n. (6) 

The equilibrium values m and n have length S cos e and 
S sin 9 respectively, and the resulting precession frequen- 
cies are 

huj+ = ^/H^ + h{js cose + h), 
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huj- = ViH - JS' sin 61)2 + h{JS cos d + h). (7) 

The slower mode involves the transverse components of n 
and the longitudinal component of m. Its frequency de- 
termines the energy gap: A = Tloj- . Figure^lshows an es- 
sentially perfect agreement between the classical analysis 
and the numerical diagonalization. This is not surpris- 
ing: because the Lieb-Mattis model has infinite-range 
interactions, the mean- field solution Q becomes exact 
in the thermodynamic limit. 

Next we explain the salient features seen in the de- 
pendence A(iJ) (Fig. O, namely the initial increase 
A (X i?^/2 and a dip around the saturation field He — JS 
for a small enough c. 

1. Weak uniform field: H <C JS 

Well below saturation, H <^ H^. = JS, the dominant 
effect is the breaking of the axial symmetry by the trans- 
verse field h. To estimate the resulting energy gap, we 
neglect the longitudinal field H and thus obtain the sec- 
ular equation 

- ?i2tj2^m = h X [(Jn + h) X (5m]. (8) 

It yields the energy gap 

A^hLo = y/h{Jn + h). (9) 

Spontaneous breaking of the axial symmetry implies 
that the staggered magnetization n attains the maximal 
length n = S for an arbitrarily weak staggered field h. 
Therefore 

A - {JShy/^ (10) 

as h — > 0. R ecalling that h — cH we obtain the gap 
A - VJScH at low fields H, in accordance with Eq. 

2. Uniform field at saturation: H — JS 

For small enough c, the gap A (if) has a minimum near 
the saturation field He = JS. Its origin can be traced to a 
reduced response of the spins to the transverse staggered 
field h at saturation. 

Inspection of the equations of motion Q in the ab- 
sence of the staggered field shows that the spin preces- 
sion can be separated into a fast mode with fko^ — JS 
and a slow mode with huj- = 0. This separation of scales 
still works in the presence of a weak transverse field. In 
contrast to the phase with a spontaneously broken ax- 
ial symmetry, the staggered magnetization n now van- 
ishes as /i — > 0. With the aid of Eq. Q, we obtain 
n = ScosO fa [2h/JSf/^S for a small h. Hence both 
terms in the expression of hu)^ in Eq. 10 are of the same 
order since H-JSsmO « [l/1){2h/ JSf/^, which yields 
the energy gap at H — JS: 

A ~ 3l/22-l/3(J^)l/3/,2/3 (11) 



Comparison of Eqs. H10|l and (|ll|l for a staggered field 
h — cH shows that the gap scales as c^/^ at weak fields 
H but becomes of order c^/^ at saturation. The latter 
is smaller (for a substantially small proportionality con- 
stant c). Therefore the initial increase of A with H will 
be followed by a dip at H ^ JS, provided that c is small 
enough (Fig. ^ . The dip disappears when c exceeds the 
critical value Cc — 0.03126(6). 



III. BEYOND THE LIEB-MATTIS MODEL 

The results presented in this paper, most importantly 
the dependence of the spin gap on the applied field H, 
were obtained for the Lieb-Mattis model. The special 
form of its Hamiltonian ^ has two important advan- 
tages. First, the existence of conserved quantities Sa 
and Sb enabled us to determine the low-energy spec- 
trum numerically for very large systems (up to iV = 2000 
spins). Second, the infinite range of spin interactions jus- 
tified the use of a mean-field approximation, allowing us 
to derive the low-energy spectrum and explain the ob- 
served features. At the same time, one must use caution 
in drawing conclusions for real magnets on the basis of 
the results obtained in an infinite-range model. In this 
section, we discuss implications of our findings for more 
realistic models — such as Eq. (QJ — in d = 1, 2, and 3 
dimensions. 

The problem of a generic Heisenberg antiferromagnet 
with two sublattices in crossed uniform and staggered 
fields does not have an exact solution. Nonetheless, the 
considerations advanced in Section 111 B\ can be extended 
to the general case of an antiferromagnet with or with- 
out long-range Neel order. To achieve this goal, we use 
a field-theoretic approach, as was done previously by 
Oshikawa and Affleck^ in their work on the Heisenberg 
chain. 



A. Effective models 

Let us identify the quantum field theory appropriate 
for the gapless state (h = 0) and in its vicinity (/i <C J). 
The uniform field H 7^ breaks the rotational sym- 
metry 0(3) down to 0(2) confining the staggered mag- 
netization n — n{cos(j),sm4),0) to the xy plane. In 
d > 2 dimensions,^** the staggered magnetization acquires 
a nonzero expectation value breaking the 0(2) symmetry. 
The low-energy degrees of freedom are long-wavelength 
fluctuations of the direction of n in the xy plane, which 
can be thought of as a Bose condensate*^ rix+iny — ne*"^. 
The spin waves are fluctuations of its phase (jj; fluctua- 
tions of the amplitude n are gapped and for this rea- 
son can be neglected. Thus one obtains an effective La- 
grangian for the low-energy excitations: 

£ = ^((/,2/s2 - |V0n -t-/OTCOS</). (12) 
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Here ps > is a spin stiffness, and s is tfie magnon ve- 
locity. The last term hn cos 0, describing the coupling to 
the staggered transverse field, breaks the residual 0(2) = 
U(l) symmetry and induces an energy gap. In d — I di- 
mension, long-wavelcngth phase fluctuations destroy the 
condensate, (ne**^) = 0; however, the low-energy theory 
(fT^ is still apphcablepi^ 

When the uniform field H reaches a critical magni- 
tude He, the antiferromagnet enters a polarized phase 
where all spins point along the field direction, as in a 
ferromagnet. For H > He, the equilibrium magnitude of 
the staggered magetization n vanishes and the effective 
field theory (|12|l no longer applies: amplitude fluctua- 
tions become soft a.t H — He- Above the critical field, 
the magnon spectrum acquires a gap S oc H — He and 
the energy dispersion switches from linear to quadratic, 
as in a ferromagnet: eic = <J + k^/2m. In the boson 
language, H = He can be viewed as the point of Bose 
condensation for magnons.^^ (There is, in fact, an exact 
mapping between spins and hard-core bosons^^ in the 
case of S" = 1/2.) The low-energy effective theory de- 
scribing universal properties of interacting bosons near 
the condensation point has been discussed by Fisher et 
alm^ It has the Lagrangian 

where ^ = Ux + iriy. When 5 <Q, the bosons condense, 
$ — ne**^, and the effective field theory reduces to that 
of phase fluctuations (fT^ . 

Next we discuss the properties of the effective mod- 
els below and at the condensation point, with particular 
emphasis on the energy gap induced by the staggered 
transverse fleld h. 



B. Weak uniform field: H He 

Well below the condensation point, phase fluctuations 
are the dominant excitations. The effective field theory 
is given by the Lagrangian shown in Eq. p2|l . Analyti- 
cal continuation to imaginary times t — —isr yields the 
classical XY model in the ordered phase in + 1 dimen- 
sions, whose properties are well known. In particular, 
the symmetry breaking field h creates a finite correlation 
length ^ ^ (h/ ps)^^^"^ ■ This translates into an energy 
gap A ^ c{h/ psy-/"^ in the quantum case. The scaling of 
the gap is the same as in the Lieb-Mattis model pO(l . 

In d = 1, the long-range order is absent. The ground 
state of the quantum model in zero transverse field corre- 
sponds to the Kosterlitz-Thouless phase with power-law 
correlations (e"'^('')e*'^^°-') ^ C /r^ with a nonuniversal ex- 
ponent 77 > 0. By the standard scaling argument, the gap 
opens as A oc /i^/^^^''/^). For the S = 1/2 Heisenberg 
chain in the limit of zero uniform field, Oshikawa and 
Affleck^ find that A cx /i^/^. 



C. Uniform field at saturation: H — He 

Near the condensation point H = He one must use the 
more complete theory (|13|l taking into account fluctua- 
tions of the condensate magnitude |$| = n. This critical 
theory has been analyzed by Fisher et al. in the con- 
text of the superfluid-insulator transition>ii Analytical 
continuation to imaginary time t = —it yields a classical 
field theory in + 1 dimensions; the singular part of the 
free energy density has a scaling form 

fsis,h) - s-^'^+'^/y'x{hs-y''/y'), (i4) 

where ys and yh are the RG eigenvalues of the "uniform 
field" S and the staggered field h; z is the dynamical 
critical exponent. 

At the saturation field He ((5 = 0), the spin gap scales 
as a power of the transverse field: A cx /i^/^'* . The expo- 
nents z and yh are most readily obtained from the trans- 
verse spin correlation function, whose long-wavelength 
part has the scaling form 

{^*ir,r)mO))- ^ (15) 

The transverse spin correlations can be computed exactly 
in the polarized phase (H > He), where the ground state 
is trivial and elementary excitations are magnons with a 
quadratic dispersion: 

(5+(r,t)5"(0,0)) cxt-'^/^exp (-^St-'-^y (16) 

where m is the magnon mass. Hence z = ys = 2 and 
yh — 2 + d/2. The S — h = fixed point is Gaussian in 
any dimension>i2. The upper critical dimension is dc = 2. 

These considerations yield the following dependence 
of the spin gap at the critical value of the uniform field 
H = He- In d = 1, A cx h^/^. In d > de = 2, the 
critical exponents revert to their mean-field values, so 
that A cx K^/^, as in the Lieb-Mattis model. 

To test the validity of the field-theoretic argument, we 
have computed numerically the energy gap A as a func- 
tion of the staggered transverse field h in & S = 1/2 
Heisenberg chain. Previous authors have examined A(/i) 
in the absence of the uniform field H^^ but not at the 
saturation point H = He — J. We have employed the 
DMRG method for system sizes up to iV = 100.^^ The 
number of states kept was set at 150 for the warm-up 
phase and 300 for the zip iteration. With the excep- 
tion of the gapless point h = 0, size effects are rather 
weak. In a chain with = 100 sites and J = 1 we find 
A — 1.60h'~'^^^ (Fig- El- The numerical results are in rea- 
sonable agreement with the predicted gap exponent of 
4/5. 

As expected on physical grounds, the energy gap A(/i) 
opens more slowly at the critical uniform field He than 
in zero field. If the staggered field is proportional to the 
uniform one, h = cH, the gap is a quantity 0{c^^^) when 
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FIG. 5: The energy gap A as a function of the staggered 
transverse field h in a S = 1/2 Ifeisenberg chain of length 
= 100 in saturating uniform field H = He = J- The 
squares are numerical results obtained by the DMRG method. 
The straight line is the best-fitting power law A — ah^ with 
& = 0.81. 

H is small and 0{c^^^) at saturation field He (in d> 2). 
Consequently, for a small enough c, the gap will have 
a local minimum near H — He, as in the Lieb-Mattis 
model. 



IV. CONCLUSIONS 

The model of a Heisenberg antiferromagnet in a uni- 
form magnetic field H perpendicular to a staggered mag- 
netic fields h — cH is relevant to real magnets with a 
staggered g tensor or staggered Dzyaloshinskii-Moriya 
interactioni^ We have determined the low-energy spec- 
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trum of such a magnet with long-range order in the model 
of Lieb and Mattis The existence of conserved quan- 
tities beyond the total spin allowed us to determine the 
spectrum numerically in large systems (up to iV = 2000 
spins). By utilizing the infinite range of interactions in 
the model, we have also reproduced the spectrum ana- 
lytically and found essentially perfect agreement with the 
numerical results. 

We have determined that the energy gap, caused by 
the breaking of the 0(2) spin- rotation symmetry by the 
induced staggered field, scales as (c77)^/^ at low uniform 
fields and is of order {cH)"^/^ at the saturation field H — 
He- The gap has a local minimum near H^ if c is small 
enough (c < 0.03126). Such a dip has been previously 
observed in numerical studies of spin chains^iS but has 
not been explained. 

Finally, we have presented scaling arguments estab- 
lishing the same power laws for the gap, and hence the 
existence of a saturation dip at small c, to any Heisen- 
berg antiferromagnet in magnetic field with long-range 
order in d > 2 dimensions. We note that a local mini- 
mum of the spin gap near the saturation field has been 
observed in SrCu2 (603)2, a Shastry-Sutherland mag- 
net with staggered g-factors and Dzyaloshinskii-Moriya 
interactions»iii 
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